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We investigate charge conductance and spin and valley polarization along with the tunnelling 
magneto-resistance (TMR) in silicene junctions composed of normal silicene and ferromagnetic sil¬ 
icene. We show distinct features of the conductances for parallel and anti-parallel spin conhgurations 
and the TMR, as the ferromagnetic—normal—ferromagnetic (FNF) junction is tuned by an external 
electric field. We analyse the behavior of the charge conductance and valley and spin polarizations 
in terms of the independent conductances of the different spins at the two valleys and the band 
structure of ferromagnetic silicene and show how the conductances are affected by the vanishing 
of the propagating states at one or the other valley. In particular, unlike in graphene, the band 
structure at the two valleys are independently affected by the spin in the ferromagnetic regions and 
lead to non-zero, and in certain parameter regimes, pure valley and spin polarizations, which can be 
tuned by the external electric field. We also investigate the oscillatory behavior of the TMR with 
respect to the strength of the barrier potential (both spin-independent and spin-dependent barriers) 
in the normal silicene region and note that in some parameter regimes, the TMR can even go from 
positive to negative values, as a function of the external electric field. 

PACS numbers: 73.23.-b, 73.63.-b, 72.80.Vp, 75.76.+j 


I. INTRODUCTION 

A close cousin of graphene, silicene has been attract¬ 
ing a lot of attention in recent times both theoretically 
and experimentally^^^, due to the possibility of new ap¬ 
plications, given its compatibility with silicon based elec¬ 
tronics. Although earlier theoretical analyses®^^ had pre¬ 
dicted the possibility of silicene and even germanene, 
stanene, etc, interest in this subject rose after experi¬ 
mentalists observed hexagonal structure in silicene sheets 
deposited on a silver substrate^^"^. Unlike graphene, sil¬ 
icene does not have a planar structure; instead it forms 
a buckled structure due to the large atomic radius of sili¬ 
con, resulting in a band gap at the Dirac point®. Further, 
it turns out that such a band gap is tunable by an ex¬ 
ternal electric field applied perpendicular to the silicene 
sheet^’®. Thus, from the point of view of applications, 
one of the drawbacks of graphene in making transistors 
is overcome by silicene, and very recently, a silicene based 
transistor has been experimentally realised®. Therefore, 
with the possibility of silicene based electronics, experi¬ 
mental interest in silicene remains quite high. 

Theoretical interest in silicene soared when it was 
realised that one could have topologically non-trivial 
phases in silicene, tuned by the external electric field^’®. 
Graphene and silicene have similar band structures and 
the low energy spectrum of both are described by the 
relativistic Dirac equation i.e., both have the Dirac cone 
band structure around the two valleys represented by the 
momenta K and K'. However, in silicene, the spin-orbit 
coupling is much larger than in graphene®’^’^®. This is 
an important difference which causes the Dirac fermions 
to become massive. Furthermore, because of the buckled 
structure in silicene, the two sub-lattices respond differ¬ 


ently to an externally applied electric field. This means 
that the Dirac mass term becomes tunable', and hence 
allows for the mass gap to be closed at some critical 
value of the electric field and then reopened. The phases 
on the two sides of the critical value where the gap is 
closed are different, with one of them being topologically 
trivial and the other being topologically non-trivial®’^’^^. 
Hence, under suitable circumstances, silicene can be a 
quantum spin Hall insulator with topologically protected 
edge states^®’^®. 

In recent years, spin based electronics or spintronics 
has become a prominent field of research both theoreti¬ 
cally and experimentally The upsurge of activity in 
this area is essentially due to the realisation that devices 
based on the spin degree of freedom are almost dissipa¬ 
tionless unlike those based on the charge degree of free¬ 
dom. The possibility of using silicene as a spintronic de¬ 
vice has been reported very recently in Ref. [15-19] due to 
its strong spin-orbit coupling. Another important quan¬ 
tity in spintronics is the tunneling magneto—resistance 
(TMR) which occurs at junctions between materials i.e., 
ferromagnet—normal metal—ferromagnet (FNF) junc¬ 
tions. The resistance of the junction is different for par¬ 
allel and anti-parallel spin configurations, and this differ¬ 
ence in resistance can be experimentally measured. 

In this paper, we study charge conductance and spin 
and valley polarizations along with the TMR in silicene 
junctions, in particular the FNF junctions, as the ex¬ 
ternal electric field is tuned through the system. We 
model our FNF setup within the scattering matrix for¬ 
malism®®. A similar setup has been utilised in graphene®^ 
to study TMR, while valley polarization has been stud¬ 
ied in a normal—ferromagnet—normal (NFN) junction in 
silicene®®. In Ref. 22, the author has investigated the con- 
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FIG. 1. (Color online) Schematic of the FNF junction in 
silicene in which two ferromagnetic patches (dark brown, dark 
grey) have been deposited on two sides of a normal silicene 
sheet (cyan, light grey) to induce ferromagnetism in it. Here 
L is the length of the normal silicene region. 

ductances and valley polarization of a NFN junctions in 
silicene. Spin and valley textures of the particle-hole ex¬ 
citations due to the addition of external fields is another 
important issue in silicene and recently, analytical and 
numerical results for the dispersion of the plasmons has 
been studied^^. However, conductances and TMR based 
on silicene FNF junction have not yet been considered in 
the literature. 

The remainder of this paper is organized as follows. 
In Sec. II, we present our model and band structure for 
the ferromagnetic silicene (FS). In Sec. Ill, we describe 
the scattering matrix for the FNF junction to compute 
the two terminal charge conductance (C?c), valley and 
spin polarizations [Vy^Vs) and the TMR. In Sec. IV, we 
present our numerical results for the conductances, spin 
and valley polarizations and TMR in the FNF set up for 
various parameter regimes. Finally in Sec. V, we present 
the summary of our numerical results followed by the 
conclusions. 


II. MODEL AND BAND STRUCTURE 

We study the FNF junction set up in silicene as shown 
in Fig. 1. Ferromagnetism is induced in silicene by the 
proximity effect when it is placed in proximity with a 
magnetic insulator, which we model by the following 
Hamiltonian given by^^ 

H = hvF{rjkxTx-kyTy) + {elE^-r](jXso)Tz + V{x)-h{x), 

( 1 ) 

where vp is the Fermi velocity of the charge carriers in 
silicene, eis the charge of the electron and ry, a correspond 
to the valley and spin indices and f corresponds to the 
sublattice (pseudospin) Pauli matrices. \so is the pa¬ 


rameter that specifies the spin-orbit coupling in silicene. 
Due to the buckled structure of silicene, the atoms in 
the two sublattices respond differently to an externally 
applied electric field E ^,'. Thus lE^ is the potential dif¬ 
ference between the two sublattices A and B due to this 
applied electric field where I is the separation between the 
two sublattices. Hence, the potential difference is a tun¬ 
able parameter and can be tuned by an external electric 
field^. Also, when the Fermi energy is close to the Dirac 
point, at the critical electric field E^ = rjaXsOy one of 
the valleys in silicene is up-spin polarized and the other 
one is down-spin polarised^. Here 77 = ±1 denotes the K 
and K' valleys respectively and tr = ±1 denotes the spin 
indices. V(x) denotes the profile for the potential barrier 
in the normal silicene region and h{x) corresponds to the 
exchange splitting, or the energy difference between the 
up and down spin electrons in the FS regions. Note how¬ 
ever, that in real materials the proximity of a ferromagnet 
to silicence can actually change the band structure of sil¬ 
icene itself. In that case, the only way to proceed will be 
to perform first principles calculations as have been done 
in graphene^^. 



FIG. 2. (Color online) Schematic of the configurations with 
P (mm) and AP {ud) spin polarizations for one of the valleys 
of the FNF silicene junction. 

We now consider the geometry shown in Fig. 1 and 
assume that the system is translationally invariant along 
the y direction. The interface between the normal and 
the FS are located at x = 0 and x = L where L the 
length of the normal silicene region sandwiched between 
the ferromagnetic patches. Here V(x) = U<d{x)Q{L — x) 
is the profile of the potential barrier modelled in the nor¬ 
mal silicene region and h{x) = h[asQ{—x) + CT's'Q{x —L)] 
denotes the exchange field or Zeeman field in the two fer¬ 
romagnetic regions with s = s' corresponding to the par¬ 
allel (P) configuration and s = —s' corresponding to the 
antiparallel (AP) spin configurations of magnetization re¬ 
spectively. We show the schematic of up and down spin 
in left region (for s = 1) and in right region for parallel 
(s' = 1 or P or uu) and antiparallel (s' = —1 or AP or 
ud) configurations in FNF junction for both E < h and 
E > h regions in Fig. 2. The orientation of the exchange 






















3 


field in the left region is kept fixed by keeping s = 1 
and then s' = 1 and s' = —1 in the right region cor¬ 
responds to the parallel (P or uu) and antiparallel (AP 
or ud) configuration respectively. Note that the E < h 
line crosses the same band in the third region for the P 
configuration, but crosses the other band for the AP con¬ 
figuration. As reported earlier^^, this feature gives rise 
to negative TMR in graphene. In silicene, also, the same 
feature is responsible for negative TMR which we discuss 
later in Sec. IV. 



k k k k 

FIG. 3. (Color online) (a) to (d) gives the schematic of the 
band structnre for s = -fl, at K (for both b and 4, spin) and 
K' (for both f and spin) valleys for ferromagnetic silicene 
for four different values of the dimensionless parameter E^/h 
(red (0.1), green (0.5), blue (1.0) and magneta (2.0)). On 
the other hand (e) to (h) gives the same for s = — 1. We 
use these diagrams to qualitatively explain the dependence of 
conductances on the electric field as mentioned in the text. 

In Fig. 3 we demonstrate the spin polarization of both 
the K and K' valleys for FS with the magnetisation di¬ 
rections up {u defined by s = -fl) and down {d defined 
by s = —1) via an energy band diagram. The diagram 
clearly shows that different spin orientations behave dif¬ 
ferently in the K and K' valleys at different values of 
the exchange field h. To visualize this picture, we fix 


^so/h = 0.5 and show the dispersion for K f, K K' f 
and K' for four different values of the tunable param¬ 
eter E^jh, for the cases u and d. Note that we use f, 
I to denote the spins of the incoming (also reflected and 
transmitted) charge carriers and we use u, d to denote the 
orientation of magnetic exchange. For the it or s = -|-1 
case, at E^/h = Xso/h = 0.5 both K f and K' have a 
vanishing band gap. On the other hand, for E^/h = 0.1, 
Ez/h = 1, and E^/h = 2, the valleys at K f, K K' f, 
K' I are all gapped. Also note that due to the exchange 
splitting h, the K' valley is shifted upwards for the | spin 
while the K valley is shifted downwards for the f spin. 
The case for d or s = — 1 is the other way around. Hence, 
it is clear that unlike in graphene, the contributions to 
the conductances from various spin configurations will 
not be identical for the K and K' valleys. 

III. THE SCATTERING MATRIX 


We model our FNF setup within the scattering ma¬ 
trix formalism^^ where we match the wave functions at 
each ferromagnet—normal interface to obtain the scat¬ 
tering matrix and find the conductances and the TMR . 
The wave functions for the valley rj in each of the three 
regions, a;<0, 0<a;<L and x > L can be written as 


^ikixX 

^ ykiSiE''^^^ 

“V2Er, 1 

^—ikixX 

1 h 1 

V R 

' —r]kiSie~ 


V A 


( 2 ) 


where oi = 1,= r for x < 0 , 02 = a, 62 = b for 
0 < X < L and 03 = t, 63 = 0 for x > L. Note that we 
keep track of the sign of the charge carriers by including 
the index Si in all the regions (the sign of the charge 
carriers changes from electron-type to hole-type when Si 
is negative in any region). This actually happens for the 
anti-parallel configuration when energy of the incident 
charge carrier is below the induced magnetic field energy 
i.e. {E < h). This charge reversal actually qualitatively 
changes the conductances, as was shown in graphene^^. 

We obtain the scattering matrix both for E > h and 
E < h hy matching the wave functions (see Eq.(2)) at 
X = 0 and x = L and solving Eq.(3) numerically. 
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Further, 

ki = \J Ei^ - {elE:^ - rjaiXsof 
and Ti = E, - {elE^ - rjaiXso) , 


with ai = a, as = a\ El = E + sh, Es = E + s'h and 
Si = sgn[Ei — {elE^ — rjaiXso)]- Since momentum is 
conserved in the y direction and does not change, it is 
convenient to write the x-component of the wave-vectors 
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as 

kxi = \Jki^ - ky^ ■ (5) 

where ky is the conserved momentum in the y direction. 

In Eqs. (4 and 5) we consider i = 1, 3 only. For sim¬ 
plicity we assume that = 0 in the middle region which 
makes the momentum in the middle region independent 
of valley and spin ((T 2 ). Hence for the central region, 

k, = ^Ei-Xl^ , 
kx2 — yj ^2 ’ 

and T 2 = E 2 ■ (6) 

where E 2 = E — U and U is the height of the potential 
barrier in the normal silicene region. 


IV. NUMERICAL RESULTS 

In this section we present our numerical results for the 
FNF junction for different parameter regimes. 

We first study the model described in Eq.(l) which 
has a spin-independent barrier in the normal silicene re¬ 
gion. We compute the conductance using the transmis¬ 
sion coefficients obtained in Eq.(3) for both the parallel 
(s = s') and the anti-parallel (s = —s') configurations 
of spins, using the Landauer—Buttiker formalism^°. We 
use the scattering matrix to compute the total transmis¬ 
sion probability T"" {di) = \t\'^kxs/k^i for parallel and 
anti-parallel configurations by choosing the spins appro¬ 
priately and for a particular incident angle 0i (which fixes 
the angles in the other regions as well). The factor of 
kxs/kxi in the transmission function is needed because 
the probablity flux density includes a factor of the veloc¬ 
ity which is essentially hk/m. Since experimentally, it is 
not easy to control the angle of incidence of the impinging 
electron, we then compute the conductance by integrat¬ 
ing over the possible angles of incidence and multiplying 
by the number of modes within the width W of the sil¬ 
icene sample. At zero temperature, this leads us to a 
conductance given by 

G"" = - -^ / T"" (6»i)cos6>id6»i . (7) 

h n Jo 

Here, Oq is the critical angle of the incident particles 
which is needed to ensure propagating particles in the 
first and third regions and is given by 9c = 7r/2 for 
ki < ko and 6c = arcsin(A:i/A: 3 ) for ki > k^. Note 
that unlike the case for graphene, in silicene, the con¬ 
tributions at the two valleys are not identical and hence, 
we do not get the degeneracy factor of two. Instead, 
the contributions at both the valleys have to be com¬ 
puted independently and added to obtain the total con¬ 
ductance through the junction. Thus we define the total 
charge conductances Gc, valley and the spin polarizations 
{Vvt'Ps) and TMR through the FNF junction in terms 



- ud 

(m) 

I 

I™ 

\ (n) 

\ 

I 

^ — uu 

GrTT _ ud 

(o) 

I 

r. — uu 

*^K-U _ ud 

(P) 

\ 

0.5 

0 

A 

0.5 

0 


0.5 

0 



E^/h EJh EJh EJh 


FIG. 4. (Color online) Conductances {Gvo-a') in units of 
e^W/Trh, for the P (uu) and AP (ud) configurations of a FNF 
junction are shown as a function of the dimensionless param¬ 
eter Ez/h for E > h [upper panels, (a-h)] and E < h [lower 
panels, (i-p)] respectively. Here Ez is the external electric field 
and h is the ferromagnetic exchange field. The value of the 
other parameters are chosen to be Xso/h — 0.5, U/h — 30. 
Energy of incident electron, for E > h is E/h = 4.0 and for 
E < h is E/h = 0.5. 


of the following constituent conductances— Gy^^,. Here, 
ss' denotes uu (P) or ud (AP) spin configurations, V de¬ 
notes the valley (K or K') and a denotes the spin of the 
incoming charge carrier in region 1 and a' denotes the 
spin of the outgoing charge carrier in region 3, which can 
be different, because we have spin-orbit coupling in the 
system. These conductances have been shown in Fig. 4 
for both E > h and E < h. It is now easy to realize, 
in reference to the band diagrams given in Fig. 3 that 
the conductances go to zero when there is a gap in the 
density of states either in region 1 or 3. The maxima can 
also be understood by noting that the density of states 
at those values of E^/h are maximum and reduce both 
when Ez/h is reduced or increased. This can be checked 
for each of the various conductances on a case by case 
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FIG. 5. (Color online) Total charge conductance (Gc) in units 
of e’^Wjnh, valley polarization (Vv) and spin polarization 
(Vs) for P and AP configurations of a FNF junction are shown 
as a function of the dimensionless parameter E^/h for E > h 
[left panels, (a), (c) and (e)] and E < h [right panels, (b),(d) 
and (f)] respectively. The value of the other parameters are 
chosen to be the same as in Fig. 4. The insets emphasize that 
Vv and Vs are actually different in magnitude for the uu and 
ud configurations for E > h regime. 


basis. 

The total charge conductance and the valley 7^®® 
and spin polarizations for both the P (s = s') and 
AP (s = —s') configurations are now defined as 


G 


ss' 

c 


V 


ss' 

V 


and Vf 


E ^ss' 

'^Vcrcr' ; 

E f/^ss' /^ss' \ 

(7(7'\^'Kg(7’ ^¥L'<7<7') 

Gf 

E f/^ss' _ /^ss' \ 

nss' 


( 8 ) 


Note that the indices ss' in Eq.(8) gives rise to four possi¬ 
ble spin conhgurations uu, ud, du, dd for a FNF junction, 
of which uu and dd imply P conhgurations with s = s' 
and ud and du denote AP conhgurations with s = —s'. 

Now from the knowledge of all the possible conduc¬ 
tances, one can dehne the tunneling magnetoresistance 
(TMR) through the FNF geometry as 


TMR = 


G““-G 

(^UU 


ud 

c 


(9) 


Note that the standard dehnition of TMR has G“‘^ in 
the denominator. However it is also sometimes dehned 
with G““^® and we choose this dehnition because in our 
case, G“®* vanishes at = h. This implies a singularity 
in the TMR which is avoided in our dehnition. Note 
that for Ez > h, the difference in TMR between the two 
dehnitions is negligible. For E^ < h, there are numerical 




FIG. 6. (Golor online) TMR is shown as a function of the 
dimensionless scale E^/h in panel (a) with green and purple 
lines corresponding to E > h and E < h regime respectively. 
In panel (b) TMR is shown as a function of U/h for E < h for 
three values of Ez/h. The other parameters are chosen to be 
the same as in Fig. 4. In the inset of panel (a) we emphasize 
the very slow rate of increase of TMR in the E < h regime. 


differences, but no qualitative difference in the behavior 
of TMR with the two dehnitions. 

The results are different for the energy regimes E > h 
and E < h, because of the difference in band structure, 
which has band gaps and hence no propagating states 
available for transport (see Fig. 3), for certain ranges of 
Ez/h for E < h. Since the main difference of silicene from 
graphene is the fact that the gap in silicene is tunable by 
the external electric held E^, we choose to focus on the 
dependence of conductances on E^. In cases, where we 
study the conductances as functions of other parameters 
such as the barrier strength U or the exchange splitting 
A, we present our results for three different values of 
Ez/h. 

Note that when E^ = Xgo, silicene is actually copla- 
nar, i.e. the two sublattices are in same plane like in 
graphene. But the spin-orbit coupling in silicene is much 
stronger than in graphene®. This increases or decreases 
the momentum of the incident charge carrier (see Fq.(4)) 
depending on the spin-polarization of the ferromagnetic 
silicene. Hence, we do not expect to reproduce the re¬ 
sults of FNF junctions in graphene in the gapless regime. 


A. Spin-independent Barrier 

Here we discuss the case where we have a finite spin- 
independent (scalar) barrier in the normal silicene region. 
The energy of the incident electron can be in the regime 
E > h 01 E < h and we present the behavior of the 
conductances and the TMR below. To carry out our 
numerical analyses, we have chosen to normalise all our 
energy scales by the Zeeman energy h, so that all our 
results are in terms of dimensionless quantities. We also 
choose to measure conductances in units of e^W/nh. 
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We first present the results for the various constituent 
conductances for the P and AP conhguration in Fig. 4 in 
order to understand the behaviour of the conductances, 
the valley and spin polarizations and the TMR in var¬ 
ious parameter regimes. The conductances are shown 
independently at the K and the K' valleys as well as in¬ 
dependently for the incoming (u) and the outgoing (a') 
spins of the charge carriers. In Figs. 4(a-h), we show the 
behaviour of the conductances at the K and K' valleys, 
with respect to the dimensionless parameter Ez/h for the 
four possibilities (tti ti) iti'I'i) in the E > h regime for 
both the uu and ud configurations. Here the uu configu¬ 
ration corresponds to the majority spin density of states 
in the left and right FS regions being up spin (parallel to 
each other) and the ud configuration corresponds to the 
majority spin being opposite (anti-parallel) in the two re¬ 
gions. When t charge carrier comes in from the left then 
it can either go to f state or j, state in right region. So 
tt etc, denote spins of the incoming and scattered charge 
carriers. The behaviour of the various conductances in 
Figs. 4(a-h) can now be understood easily when analysed 
in terms of the band structure presented in Fig. 3. 

For E > h, the results have been presented for E/h = 
4. Using the band diagram, it is easy to check that both 
at K and K', there are always electron states available 
for conductance for both the P and AP configurations 
and for all possible incoming and outgoing spins. The 
differences in the magnitude both at K and K' valleys 
stems from the decrease in the momentum of the prop¬ 
agating states at the Fermi energy, as can be seen from 
the band diagram. 

The results ior E < h is shown in Figs. 4(i-p) for 
E/h = 0.5. Consider the uu case for the valleys K 
and K'. For GK.K'tt? shown in Fig. 4(i) and 

Fig. 4(m)), the band diagram (see Fig. 3) shows that if 
we start with a spin up electron at the K valley, (Fig. 4(i)) 
then there is a non-zero density of states for f electrons 
in the third region for all values of E^/h until it reaches 
the value of 2 (the magenta line goes above the value 
oi E/h = 0.5). On the other hand, for the K' valley, 
(shown in Fig. 4(m)), there is no density of states for 
the t electrons beyond E^/h = 1.0 (the blue line goes 
above the value E/h = 0.5) in the third region. This 
explains why beyond E^/h = 2 for the K valley and be¬ 
yond Ez/h = 1 for the K' valley, the conductances 
and are zero. It is also clear from the band dia¬ 

gram that for G^.^, its value increases from the value 
at Ez/h = 0, because the momenta of the electrons at 
E/h = 0.5 grows (comparing the red and green lines) and 
beyond that it decreases (comparing the green, blue and 
magenta lines). On the other hand, for it is clear 

that the momentum of the electrons at E/h = 0.5 de¬ 
creases as a function of Ez/h (comparing the red, green 
and blue lines). This explains why the conductance rises 
initially and then falls beyond Ez/h = 0.5 for gJ^ and 
why it falls monotonically for G^,. 

A similar detailed analysis can also be made for the ud 
case as well as each of the other graphs in Figs. 4(j,k,l) 


and Figs. 4(n,o,p), which explains each feature of the 
graph. However, since the method is similar to what 
has been described above, we will not go through each 
one of the graphs in detail. The behaviour of the charge 
conductance, the valley and spin polarizations and the 
TMR are also now understandable, since we can explain 
how each of the constituents Gy^^, behave as a function 
of Ez/h from the band diagram. 
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FIG. 7. (Color online) In the left panels, we show the total 
charge conductance (Gc) in units of e^W/irh, valley and spin 
polarizations ('P^'* , Vs" ) for E > h and on the right panels 
for E < h respectively. The dashed and solid lines correspond 
to two different exchange splittings (A//i) due to the spin- 
dependent barrier in the normal silicene region. The value of 
the other parameters are chosen to be Xso/h = 0.5, U/h = 
30. Energy of the incident electron, for E > h, is E/h = 4.0 
and for E < h, E/h = 0.5. The small difference in magnitude 
of Vv" and Vs" is highlighted in the insets for E > h regime. 


{a) E > h 

In Figs. 5(a), 5(c) and 5(e), we show the behaviour of 
the charge conductance (G“ ) and the valley and spin 
polarizations (V/" and Vg" ) with respect to the dimen¬ 
sionless parameter Ez/h for the P and AP configurations 
{uu and ud) in the E > h regime. Note that in Fig. 5(a), 
G““ and G/^’’ are both finite at Ez/h — 0 and start de¬ 
creasing as we increase the value oi Ez/h. Gc is obtained 
by summing G^ and G||/, which in turn are obtained 
as 

G ss' f/^ss' I ^ss' \ I (/^ss' I ^SS N 

K — V^Ktt 

+ ( 10 ) 

g' 

and similarly for G^,. In other words, the total charge 
conductance G"" = G^ -I- G^, is obtained by summing 
over all the conductances in the panels (a) to (h) in Fig. 4. 

In Fig. 5(c), "P"" and are plotted which are close 
to zero on the scale of the charge conductance. How- 
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ever, they are not identical, as shown in the inset. But 
it appears that in this regime, silicene has negligible val¬ 
ley polarization, similar to graphene, which in fact has 
no valley polarization at all, since the two valleys are 
identical. This can be understood because the valley po¬ 
larization Vv is simply proportional to G|| — Gf|,, and 
as can be seen from Fig. 4 that the magnitudes of the con¬ 
ductances at K and K' are almost the same for E > h. 
In Fig. 5(e), the behaviour of the spin polarization has 
been shown for both P [uu) and AP [ud) configurations, 
which is also very small in this regime. As shown in the 
inset, the spin polarization is positive for the P and nega¬ 
tive for the AP configurations and increases as a function 
of Ez/h. This difference is due to the spin-orbit coupling 
in silicene, whereas in graphene, they are much smaller, 
since the spin-orbit coupling is vanishingly small. Fi¬ 
nally, in Fig. 6(a), the behavior of the TMR is shown as 
a function of E^/h by the green solid line in the E > h 
regime. In this regime, the TMR is close to zero. 

(h) E <h 

The right panels in Fig. 5 shows the behaviour of the 
charge conductance G®^ and the valley and spin polar¬ 
izations Vy^ , with respect to the dimensionless pa¬ 
rameter Ez/h^ in the E < h regime, for the same spin 
and polarization configurations mentioned earlier. These 
are just the appropriate sums and differences of the con¬ 
stituent conductances in Figs. 4(i-p). Here also, their be¬ 
haviour is easy to understand by comparing each of the 
graphs in Figs. 4(i-p) with the band diagrams in Fig. 3 
and noting when there is no density of states for the con¬ 
figuration in either the incoming or the outgoing spin 
configuration of the charge carriers. For instance, for the 
uu case, there is only one contribution for Ez/h > 1 in 
Fig. 4(i) and for the ud case, there is no contribution 
for Ez/h > 1. This is because we have chosen E/h = 0.5 
and the blue line {Ez/h = 1.0) in the band diagram goes 
above that line either for the incoming or scattered re¬ 
gion for all cases in the ud configuration and all but one 
case in the uu configuration. In other words, their be¬ 
haviour follows what is expected from the availabality or 
non-availability of propagating states at the K and K' 
valleys as explained above in the discussion of Fig. 4. 

The most interesting point to note is that the valley 
polarization and the spin polarization for the parallel or 
uu configuration is unity for Ez/h > 1 in E < h regime. 
This is simply because of the entire contribution to the 
conductance in this regime originates from G^^^. So 
the conductance is both fully valley and spin polarised 
and would be an important regime to achieve by tuning 
the incident electron energy E/h < 1 and the electric 
field Ez/h > 1. In the anti-parallel or ud regime, the 
spin polarization can be tuned to negative values when 
Ez/h < 1, but without any valley polarization. 

The behaviour of TMR is demonstrated in Figs. 6(a) 
and 6(b) for the E < h regime. For Ez/h = 0, the TMR 
is negative and reaches its maximum negative value. 
Then the TMR increases as we increase the value of Ez/h 



E^/h U/h U/h 

FIG. 8. (Color online) The variation of the TMR with respect 
to Ez/h and U/h is shown for two values of A/h. In the 
panels, the dashed and the solid lines correspond to A/h = 0.5 
and A/h = —0.5 respectively. We choose the same value of 
the other parameters as mentioned in Fig. 7. 


and changes sign and reaches saturation at Ez/h ~ 1 
since the conductance becomes fully spin polarised at 
that point. Note that we need to restrict the value of 
Ez/h below two as the TMR takes an indeterminate form 
at Ez/h = 2 due to the vanishing of both G““ and G“'^ 
(see Fig. 5(c) and Eq.(9)) for all spin configurations. This 
is a consequence of our choice of the incident energy at 
E/h = 0.5. 

The striking feature of positive to negative transition 
in the TMR also arises as we vary the strength of the 
potential barrier U in the middle normal silicene region 
for different values of Ez/h. This feature is shown in 
Fig. 6(b) where the TMR oscillates between positive and 
negative values with respect to U for Ez/h = 0.1. Such 
oscillations of the TMR from positive to negative values 
have been reported earlier in Ref. 21 in graphene, due 
to the change in the type of the charge carrier in the 
third region. Note also there is no significant qualitative 
change in the behavior of TMR even if we choose U ~ 
h. This extra tunability of the TMR with respect to an 
external electric field is a unique feature of silicene that 
we wish to emphasize here in this manuscript. 


B. Spin dependent Barrier 

Here we discuss the effect of a spin-dependent barrier 
on the total charge conductance, valley and spin polar¬ 
izations and the TMR. The barrier is modelled in the 
normal silicene region as Ua- = U — a A which is shown 
in Fig. 2. Here positive (negative) A represents the ex¬ 
change splitting in the silicene barrier with its magneti¬ 
zation parallel (anti parallel) to the spin orientation of 
the FS in the first region. 

In Figs. 7(a-f) we show the behaviour of the total 
charge conductances G®® , spin and valley polarizations 
7^1® and 7^^® , in the E > h and E < h regimes, for 
±A//i. Since the qualitative behavior of all the con- 













ductances remain similar to the spin independent bar¬ 
rier case, we do not show the behaviour of the conduc¬ 
tances at the K and K' valleys independently, or analyse 
the graphs in detail via the band structure. Similarly, 
in Figs 8(a-c), we show the behaviour of the TMR as a 
function of Ezjh and as a function of U/h as well (for 
Ez < h), for both ±A//i. We find that the results are 
fairly similar to the spin-independent barrier case. 

V. SUMMARY AND CONCLUSIONS 

To summarize, in this paper, we have investigated the 
transport properties (charge conductance as well as spin 
and valley polarizations) and the TMR through a FNF 
junction in silicene. Here we have adopted the Landauer- 
Buttiker formalism to carry out our analysis. We show 
that the conductances and the TMR in this geometry can 
be tuned by an external electric field E^ for each case (ftj 
ti, it; it) in the left and right ferromagnetic silicene re¬ 
gions, for both parallel (uu or dd) and anti-parallel {ud 
or du) configurations. For specific values of the electric 
field, we analyse both the charge conductance and valley 
and spin polarizations in terms of the independent be¬ 
haviour of the conductances at the two valleys and the 
band structure at specific incident energies. We find that 
we can tune a fully valley polarised and also a fully spin 
polarised current through our setup via the external elec¬ 
tric field. We also find that the TMR can be tuned to 100 
% in this geometry via the electric field. This is one of the 
main conclusions of our analysis. We also show that the 
TMR through our setup exhibits an oscillatory behavior 


as a function of the strength of the barrier (both spin 
independent and spin dependent) in the normal silicene 
region. The TMR also changes sign between positive and 
negative values and such a transition can be tuned by the 
external electric field. This is another conclusion of our 
analysis. Hence, from the application point of view, our 
FNF geometry may be a possible candidate for making 
future generation spintronic devices out of silicene. 

As far as the practical realization of such a FNF struc¬ 
ture in silicene is concerned, it should be possible to fabri¬ 
cate such a geometry with the currently available exper¬ 
imental techniques. Ferromagnetic exchange in silicene 
may be achieved via proximity effect using a magnetic 
insulator, for instance, EuO^'*’^®. The typical spin or¬ 
bit energy in silicene is Xso ~ 4 meV®. For an incident 
electron with energy A ~ 4 meV and exchange energy 

~ 8 meV, the maximum value of the spin and valley po¬ 
larization as well as sign change of TMR from positive to 
negative value occur at an electric field E^ ^ 0.03 VA“^, 
potential barrier of height U ^ 160 meV, exchange split¬ 
ting in the normal silicene region A ~ 4 meV and width 
of the barrier L ^ 100 nm. 
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